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d ' In the framework of RGM, the binding energy of one channel AA(3o)(d*) and AA(o,3) 

are studied in the chiral SU{3) quark cluster model. It is shown that the binding energies 
of the systems are a few tens of MeV. The behavior of the chiral field is also investigated by 
comparing the results with those in the SU{2) and the extended SU{2) chiral quark models. 
It is found that the symmetry property of the AA system makes the contribution of the 
relative kinetic energy operator between two clusters attractive. This is very beneficial for 
forming the bound dibaryon. Meanwhile the chiral-quark field coupling also plays a very 
important role on binding. The S-wave phase shifts and the corresponding scattering lengths 
of the systems are also given. 



Abstract 



^This work was supported in part by the National Natural Science Foundation of China under grant number 
1977551B and the Chinese Academy of Sciences under grant number B78 



1 Introduction 

With the consideration that the Quantum Chromodynamics (QCD) is the underlying theory 
of the strong interaction, the hadronic systems are assumed to be composed of quarks and 
gluons, and most of descriptions of baryons are associated with some effective degrees of 
freedom of nonperturbative QCD (NPQCD), such as constituent quarks [1], and chiral fields 
[2], etc.. But, in the multi-baryon system, the property of the system mostly can be well 
explained in the baryon and meson degrees of freedom with only some special cases where the 
quark and gluon degrees of freedom has to be taken into account as exceptions [3]. In principle, 
nothing special would prohibit a system with the baryon number larger than one, and there 
might exist dibaryons which are composed of six quarks. Many efforts have been made in 
both theoretical and experimental investigations to search such object in past years [4-7]. 
Unlike the deuteron, the dibaryon has a QCD motivated origin and should be a system with 
quarks and gluons confining in a rather smaller volume, say smaller than 0.85fm in radius. 
In this system, both the perturbative QCD (PQCD) and the NPQCD effects should be taken 
into account. However, because of the complexity of NPQCD, one has to use effective models 
to simulate the strong interaction, especially in the NPQCD region. Thus, investigating 
dibaryons is a prospective field to test various models of NPQCD, and consequently, to 
enrich our knowledge of the QCD phenomenology. 

According to Jaffe's calculation with the bag model [4], the color magnetic interaction 
(CMI) of the one-gluon-exchange (OGE) potential of the H particle is attractive. In 1987, 
K.Yazaki [8] investigated non-strange six-quark systems by considering OGE and confinement 
potentials and showed that CMI in the AA(3 g) system demonstrate the attractive feature 
among the NN ^ NA and AA systems. Therefore, this state Is highly possible a dibaryon. 
Studying this system is particularly plausible. 

Since Jaffe predicted the H dibaryon in 1977 [4], many works have been done on searching 
the possible existences of dibaryons such as strangenessless dibaryons d* [5,8-10] and d' [11] 
and the strange dibaryons such as H, QQ, 'B*0, and HO [13,14,12]. The current calculations of 
H showed that the mass of H is around the threshold of the AA channel [14,15]. This result is 
consistent with the reports of current experiment [16]. The recent theoretical result of d* is 
obtained by considering the coupling of the A A and CC (hidden color) channels in the chiral 
quark models. It is shown that the binding energy of d* is in the order of tens MeV [10]. 
The predicted nn, E*n and EQ are about 100 MeV, 80 MeV and 30 MeV, respectively [12, 
23-25]. 



In the theoretical study of dibaryon, the quark cluster models have been Intensively em- 
ployed. With these models, one can easily realize the antisymmetrization character between 
quarks belonging to different clusters, which can only be neglected when clusters are well 
separated from each other. 

For the system with a radius not larger than 0.85fm, the Pauli principle on the quark 
level would play a substantial role [20]. According to our calculation, AA(3o) (d*) has 
the most antisymmetized property in the spin-flavor-color space, which is characterized by 
the expectation value of the antisymmetrizer in the spin-flavor-color space. This property 
appears in another non-strange six-quark system, AArQ^^y In this paper, we would discuss 
the one channel AA(o,3) together with d* within the chiral quark cluster model. The paper 
is organized in such a way that the model is briefly presented in Sect. 2, the results and 
discussions are given in Sect. 3 and the conclusions is drawn in Sect. 4. 

2 The effective Hamiltonian and two-cluster wavefunction 

As an effective theory of QCD, in the chiral quark model, the constituent quark and the 
chiral field are assumed as the effective degrees of freedom, and a long-range confinement 
term, which dominants the long-range NPQCD effect and still cannot strictly be derived 
up to now, a short-range one-gluon-exchange term, which describes the short-range QCD 
effect and a set of terms induced by the coupling between the quark and chiral fields, which 
mainly depict the short- and medium-range NPQCD effects are considered as the effective 
interaction. Thus the dynamics of a six-quark system in the SU(3) chiral quark model is 
governed by the effective Hamiltonian [18,19] 

H = T + E(^r^^ + V;^^^ + l/.f + V^), (1) 

i<j 

where T denotes the kinetic energy operator of the system and V^j , V^j , V^ and Vfj 
represent the confinement, one-gluon exchange, pseudo-scalar chiral field induced and scalar 
chiral field induced potentials intervening between the i— th and j— th quarks, respectively. 
The confinement potential can phenomenologically take a quadratic form 



where the superscript a is the color index. The OGE potential can be derived from the 
perturbative tree diagram 
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with the tensor operator Sij being 



3(i?i • hj) {Sj ■ hj) - i^i ■ dj] 



(4) 



The pseudo-scalar and scalar field induced potentials are originated from the restoration of 
the important symmetry of strong interaction — the chiral symmetry and can be written as 



V, 



PS 



«j 



C{gch, m^ai ^) 



m; 



12m jm 



'3 



[/i(m^^, A, Vij) {ai ■ Uj) + f2{m^a-: K nj) Sij ](A^A|)/ , 



(5) 



and 






-C{gch, in^a^ ^) • [fsim^a^ ^> ^y) 



+ 



m 



AmiTJij 



^ U{m,^, A, n,) L ■ {a- + a,) ]{X^X'^)f 



(6) 



respectively. In above equations, the subscript / denotes the flavor index. The functions 
fi, Y, G, and H and the constant c- number C are 
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respectively, with 
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In Eq.(5), tTq with (a = 1,2,. ..,8, and 0) correspond to the pseudoscalar fields vr, K, rjs 
and 770 ) respectively, and rj and ry' are the linear combinations of r/o and r/g with the mixing 
angle 6 [14]. In Eq.(6), aa with (a = 1,2, ...,8, and 0) denotes the scalar fields a', k, e and 
a, respectively. 

In the framework of Resonating Group Method (RGM), the wavefunction of the Six-quark 
system can be written as 



*6, = A[<^A'i>BXi^AB)Z(Rcn.)] (17) 

where (J)a(b) denotes the wavefunction of cluster A(B), xi^As) is the trial wavefunction 
between clusters A and B, Z(Rcm) is the wavefunction of the center of mass motion (CM) 
of the six quark system and A represents the antisymmetrizer. Expanding the unknown 
wavefunction xi^As) by well-defined basis functions, such as Gaussian functions, one can 
solve RGM bound state equation to obtain eigenvalues and corresponding wavef unctions, 
simultaneously. The details of solving RGM bound-state problem can be found in Refs. 
[21,22]. 

The model parameters should be fixed at the very beginning by the mass splittings among 
A^,A, S and S, respectively, and the stability conditions of octet (S = 1/2) and decuplet 
{S = 3/2) baryons, respectively. The resultant values are listed in Table 1. 

Table 1 Model parameters under SU{3) chiral quark model 





Sctl 


Set2 




Sctl 


Sct2 


ruu {MeV) 


313 


313 








rus [MeV] 


470 


470 








K ifm) 


0.505 


0.505 








m^ if-m.-^) 


0.7 


0.7 


Ax ifm ^) 


4.2 


4.2 


ruk {fm~ ) 


2.51 


2.51 


Ak (/m-i) 


4.2 


4.2 


rrin (/m"^) 


2.78 


2.78 


A^ (/m 1) 


5.0 


5.0 


nirj' (/m"^) 


4.85 


4.85 


A^/ (/m-i) 


5.0 


5.0 


rria ifm~^) 


3.17 


3.17 


A. (/m-i) 


4.2 


7.0 


rricr' {fm~^) 


4.85 


4.85 


A.' (/m-i) 


5.0 


5.0 


rrif, (/m-^) 


4.85 


7.09 


Ak (/m-i) 


5.0 


7.61 


me (/to-1) 


4.85 


7.09 


A. (/m-i) 


5.0 


7.61 


Qu 


0.936 


0.936 








9s 


0.924 


0.781 








Quu {MeV/Jni^) 


54.34 


57.71 


at (MeV) 


-47.69 


-48.89 


Qus {MeV/fm^) 


65.75 


66.51 


at (MeV) 


-41.73 


-50.57 


ass (MeV/fm^) 


102.97 


115.39 


at (MeV) 


-45.04 


-68.11 



In this table, rriA denotes the mass of particle A which can be either valence quark 
or meson involved, A^ represents the corresponding cut-off mass of the particle, Qq is the 
coupling constant of gluon to the valence quark q, flgigj depicts the confinement strength 
between valence quarks qi and q2, respectively, and Qg^gj denotes the corresponding zero-point 
energy. It should be mentioned again that with both sets of parameters, one can reasonably 
reproduce the NN and NY scattering data and part of single baryon properties and renders 
a mass of the H particle which agrees with the available experimental data [17-19,23-25]. 



3 Results and discussions 

We are now ready to study the structure of AA dibaryons. First of all, let us look at the 
symmetry character of the AA system. As is emphasized above, the antisymmetrization of 
quarks belonging to different clusters plays a substantial role in the structure of the dibaryon. 
In both AA(3 0) and AAm^s) systems, the expectation value of he antisymmetrizer in the 
spin-flavor-color space, {A^-''^), is equal to 2. As a result, the kinetic energy term provides an 
attractive effect. Therefore, although the CMI of AA(o,3) does not have an attractive nature 
like AA(3o), it is still possible to form bound AA(o,3) like AA(3o). 

With parameter Sets 1 and 2, we calculate the binding energy and the root-mean-square- 
radii (RMS) of two systems. The results are tabulated in Table 2. 

It is shown that both AA(3 q) and AA(o^3) are bound states, their binding energies are 
a few tens of MeV and the correspondent RMS' are around 1 fm. The binding nature of 
AA(o,3) is consistent with our speculation. The binding energy of AA(o,3) is lower than that 



of AA(3 0) by 5-6 MeV. This is attributed to the fact that the CMI of AA(o,3) is repulsive in 
contrast to that of AA(3 q), although their antisymmetrization characters are the same. The 
results with Sets 1 and 2 have small difference mainly because of the variation of Qq brought 
by the mass changes of k and e, although the strange meson clouds are not important in the 
non-strange system. 

Then, we try to reveal the effects of chiral-quark interactions in AA systems. This can 
be realized by locating calculations in other two chiral quark models, the extended SU{2) 
and SU{2) chiral quark models. In the former one (hereafter call Model II), the chiral-quark 
interactions are provided by scalar meson a and all pseudoscalar mesons vr, K, rj and rj/, 
while in the latter one (hereafter call Model III), the chiral-quark interactions are induced by 
the pseudoscalar meson vr and scalar meson a only. The results in these two models are also 
tabulated in Table 2. It is seen that the binding natures of the two systems remain the same, 
but the values of binding energies and RMS' have substantial differences. In the higher spin 
system like AA(3o), the binding energy would increase when the model changes from I to II 
or III, and the largest value occurs in Model II. But in the lower spin state like AA(q3), the 
binding energy changes in the opposite direction and the largest value appears in Model I. 
Anyway, the fact that the binding energy of the system with spin S* = 3 is larger than that 
of the system with spin S = remains un-changed when the model used is alternated. 
Table 2 Binding energy Bab and RMS for AA(3o) and AA(o 3) systems^ 



Channel 


one channel 

AA(3,o) 


one channel 

AA(o,3) 


^AA(3.o) 

(MeV) 


RMS 
ifm) 


^AA(u.3) 
(MeV) 


RMS 
ifm) 


Model I Setl 


22.2 


1.01 


16.0 


1.10 


Model I Set2 


18.5 


1.05 


13.5 


1.14 


Model II 


64.8 


0.84 


6.3 


1.25 


Model III 


62.7 


0.86 


13.2 


1.11 



f Bab denotes the binding energy between A and B baryons and RMS represents the corresponding root- 
mean-square radius. 

This result can be understood in the following way. The basic observation is that their 
symmetry structures, which cause (A)^'''^ = 2, where sfc denotes that the operation takes 
place in the spin-flavor-color space only, is very beneficial in forming bound state. Namely, 
the quark exchange effect due to the Pauli principle is enormous so that the quarks in two 
clusters can be sufficiently close, and consequently the contribution from the kinetic energy 
term shows the attractive character. On the other hand, the a field induced interaction 



provides a fairly strong attraction which plays a dominant role in forming bound state. 
These two factors make the binding nature of the system stable with respect to models and 
parameters. 

Moreover, from the form of Gell-Mann matrices we know that the K and k meson ex- 
changes do not exist between u or d quark pairs, namely the chiral fields K and k do not 
contribute in both systems. The contributions from pseudoscalar mesons in two systems are 
comparable. But the contribution from the at meson are characteristically different in these 
two systems. Opposite to the contributions of CMI, it is strongly repulsive in AA(3 g) and 
relatively weakly attractive in AA(o,3), and the repulsive strength of ai is much larger than 
the attractive strength in CMI. As a result, when one moves from Model I to Model II or 
III, the scalar mesons o"/, n and e are turned off, the wavefunction of AA(3 q) would move 
inward and that of AA(q 3) would distribute more outward. Consequently, the contribution 
from a becomes stronger in AA(3o) and a little weak in AA(o^3). That is why when model 
alternates from I to II or III, AA(3 g) becomes much more bounder and AA(o,3) turns out to 
be a little less bound. 

Above results can be crosschecked by their S-wave scattering phase shifts and correspond- 
ing scattering lengths. We plot the S-wave phase shifts of these two systems in Figs. 1 and 2, 
respectively. The solid and dashed curves are those with parameter Sets 1 and 2, accordingly. 
With these phase shifts, one can easily extract the corresponding scattering lengths, which 
are tabulated in Table 4. 

Table 4 The Scattering length a of the AA systems 





one channel 

AA(3,o) 


one channel 

AA(o,3) 


Model I Setl 


-1.30 (/m) 


-2.56 ifm) 


Model I Set2 


-1.16 (/m) 


-2.72 ifm) 



Both phase shifts and scattering lengths are consistent with our above results. 



4 Conclusions 

In the AA system, there exist two bound states, AA(3o) and [d*) AAm^s). The former one 
has been predicted as a bound dibaryon [5,8-10]. Its binding energy varies model by model. 
By employing the SU{3) chiral quark model with which the available empirical data can be 
well-reproduced, the binding energy of the one channel AA(3 g) (or d*) is 18.5-22.2 MeV and 
the corresponding RMS is about 1.0 fm. This result is consistent with most reports of other 



theoretical predictions [5,8-10]. Then, we predict the binding energy of AAj-g^s). The result 
turns out to be 13.5-16.0 MeV and the corresponding RMS is around 1.1 fm. The binding 
behaviors of these systems are attributed to their special symmetry properties which make 
the contribution of kinetic energy term beneficial for forming dibaryons. Moreover, the chiral 
fields also provide substantial contributions to their binding behaviors. Among these chiral 
fields, the relatively stronger a induced interaction is always the dominant one, no matter 
which model and model parameters are used. As a consequence, the binding phenomena of 
the AA(3 0) and AA(o,3) systems are always remained the same. 

However, the mass of these two systems predicted in the one-channel approximation are 
all larger than the thresholds of strong decay channels AA -^ ANir and AA — > NNtttt which 
are about 155 MeV and 310 MeV ( even after considering the CC channel coupling in d* , 
this conclusion is still true). These dibaryons would have very broad widths. Because the 
newly predicted AA(o,3) has the lower spin S = with respect to AA(3 g), it might be more 
favorable to be experimentally detected. 
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